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7 Abstract: We have recently reported on a theoretical digital signal-processing algorithm for improved energy and 

8 position resolution in position-sensitive, transition-edge sensor (PoST) X-ray detectors [Smith et al., Nucl. Instr and 

9 Meth. A 556 (2006) 237], PoST’s consists of one or more transition-edge sensors (TES’s) on a large continuous or 

10 pixellated X-ray absorber and are under development as an alternative to arrays of single pixel TES’s. PoST’s 

1 1 provide a means to increase the field-of-view for the fewest number of read-out channels. In this contribution we 

12 extend the theoretical correlated energy position optimal filter (CEPOF) algorithm (originally developed for 2-TES 

13 continuous absorber PoST’s) to investigate the practical implementation on multi-pixel single TES PoST’s or 

14 Hydras. We use numerically simulated data for a nine absorber device, which includes realistic detector noise, to 

15 demonstrate an iterative scheme that enables convergence on the correct photon absorption position and energy 

16 without any a priori assumptions. The position sensitivity of the CEPOF implemented on simulated data agrees very 

17 well with the theoretically predicted resolution. We discuss practical issues such as the impact of random arrival 

1 8 phase of the measured data on the performance of the CEPOF. The CEPOF algorithm demonstrates that full-width- 

19 at-half-maximum energy resolution of < 8 eV coupled with position-sensitivity down to a few 100 eV should be 

20 achievable for a fully optimized device. 

21 

22 Keywords : Transition-edge sensor; position-sensitive X-ray detector; microcalorimeter; International X-ray 

23 Observatory; correlated energy position optimal filter. 

24 PACS: 07.85. Fv; 95.55. Ka; 95.75. Pq 


* Corresponding Author, email: stephen.j.smith@nasa.gov; tel.: 1-301 -286-3719; fax: 1-301-286-1684. 


1 


1 1. Introduction 

2 At NASA Goddard Space Flight Center, we are developing microcalorimeters such as transition-edge sensors 

3 (TES’s) [1] and metallic magnetic calorimeters (MMC’s) [2] for X-ray astronomy missions such as the International 

4 X-ray Observatory (IXO), jointly under study by NASA and ESA. IXO is a merger of NASA’s Constellation-X [3] 

5 and ESA’s XEUS [4] mission concepts. Beyond IXO, this development is targeted at NASA’s Generation-X ‘vision 

6 mission’ [5]. Arrays of single pixel TES’s offer impressive full-width-at-half-maximum (FWHM) energy resolution 

7 of A£ ~ 2-3 eV in the soft X-ray energy range coupled with count-rates that can accommodate up to 100’s per pixel 

8 per second [6-8]. Alternatives to single pixel arrays are position-sensitive TES’s (PoST’s) [9]. PoST’s consist of 

9 either a single large absorber or a series of segmented absorbers, which are thermally coupled to one or more read- 

10 out sensor [10-13]. Heat diffusion between the absorber(s) and the TES’s gives rise to position dependence of the 

1 1 pulse shapes and enables position discrimination. Such devices are aimed at reducing the number of read-out 

12 channels for a given focal plane area at some expense in energy resolution and count-rate capability. 

13 

14 We are currently pursuing single-channel PoST’s, commonly referred to as Hydras, which consist of multiple 

15 absorber elements each with a different thermal coupling to a single TES readout channel. Results from first 

16 generation 4-pixel Hydras have demonstrated FWHM energy resolution of between 5-6 eV at an energy of 5.9 keV 

17 for devices with a 250 pm pixel pitch [13-15]. We are extending this approach to incorporate more absorber 

18 elements, enabling further increases in the field-of-view for missions such as IXO, and further, to investigate 

19 alternative MMC Hydra designs. 

20 

21 The ability to differentiate between different pulse shapes scales with the photon energy £ [16]. For a segmented 

22 absorber PoST, such as the Hydra, the pulse shapes are quantized and the position resolution (Ax) is constrained by 

23 the pixel size. Below some threshold energy the signal-to-noise will be degraded such that it is not possible to 

24 distinguish between different shapes of pulse and position sensitivity will be lost. For a continuous absorber device 

25 it follows that Ax cc ME. Thus continuous PoST’s can potentially achieve very high position resolution. However in 

26 the continuous devices, the position x and energy E are intrinsically correlated such that an error in thex (due to 

27 noise) will propagate into an error in E. Depending upon the specific device design, this can result in a significant 
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1 degradation in AE. The signal processing analysis is typically more computationally intense for the continuous 

2 absorber case since there is a continuous population of different pulse shapes, as opposed to a discrete set 

3 corresponding to each element of a segmented absorber. The energy of the photon is found using a digital optimal 

4 filter [17], optimised for a specific pulse shape. Consequently a much larger number of template filters will be 

5 required (of the order Lt Ax, where L is the linear dimension of the detector), otherwise the filters will not be well 

6 matched to the pulse shapes and some loss in resolution will result. 

7 

8 We are developing complex signal-processing algorithms for position-sensitive microcalorimeters [15,16,18]. In ref. 

9 [16] we reported a theoretical signal-processing algorithm for simultaneous determination of the energy and position 

10 in one-dimensional (continuous absorber) 2-TES PoST’s. In this contribution we adapt this correlated energy- 

1 1 position optimal filter (CEPOF) algorithm to investigate segmented absorber devices such as the Hydra, which we 

12 are currently developing for 1X0. We use numerically simulated pulses (for a 9-pixel Hydra), which include realistic 

13 detector noise to demonstrate an iterative implementation method, which converges on the correct photon position 

14 and energy without any a priori assumptions. 

15 

16 2. Hydra numerical model 

17 To simulate realistic pulse and noise data we use the fourth order Runge-Kutta [19] technique to integrate a system 

1 8 of ordinary non-linear differential equations that describes the Hydra. For this study we assume a 9-pixel Hydra, 

19 which is an alternative to the 4-pixel devices already experimentally characterised. Fig. 1 shows a schematic 

20 diagram of the 9-pixel Hydra concept. Each Hydra absorber is modelled as a single heat capacity C coupled to a 

21 TES via a link conductance G L . We use typical design parameters from the experimental characterisation of the 4- 

22 pixel designs. Table 1 lists the design parameters used in this study. We assume a linear current and temperature 

23 dependent resistive transition of the form 

24 

25 R{T,l) = R 0 + a^AT + P^AI, (i) 

26 
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1 where R 0 , T 0 and 7 0 are the detector resistance, temperature and current at the quiescent operating point and a and p 

2 are the unit- less TES resistance and current sensitivities respectively. Over typical resistance excursions due to 5.9 

3 keV photons, such a transition shape is found to describe the experimentally measured pulse shapes accurately. To 

4 incorporate realistic detector noise we use a random number generator with Gaussian distribution, which is added to 

5 the input of the model at every time step. For the TES near equilibrium, non-linear Johnson noise [20] we set the 

6 variance of our Gaussian noise to (4k b T 0 RoAB( 1+2/J)) 1/2 (in units of volts) where k h is the Boltzmann constant and 

7 A B is the bandwidth of the noise, which in the context of our numerical model is equivalent to 1/(2A/), where At is 

8 the integration step size. We also include the Johnson from the TES bias circuit, which has the same form except it 

9 is in equilibrium hence we drop the j3 dependence. For the thermal fluctuation (phonon) noise between the different 

10 detector heat capacities as well as the detector and the heat sink, we have (4khT 0 2 GyAB) V2 (in units of watts), where 

11 G is the thermal conductance between the relevant detector elements and yis a unit-less constant which describes the 

12 temperature gradient across the conductance link [21]. It is ~ 0.5 for the thermal fluctuation noise between the TES 

13 and heat sink and is ~ 1 for the internal thermal links which are approximately isothermal. In our 4-pixel Hydras we 

14 observed an excess in the low frequency noise of approximately 1 .2 times the level of the phonon noise to the heat 

1 5 sink [14], which we also include in this model. A white readout noise equivalent to 5 pA/VHz on the TES input is 

16 also assumed. This model intrinsically includes non-linearity of the pulse shape as well as any associated non- 

17 stationary noise. Fig. 2a shows a plot of the simulated pulse shapes for the 9-pixel Hydra with design parameters as 

18 listed in Table 1, the noise terms are set to zero. The internal thermal conductances are chosen to give approximately 

19 equal spacing between the pulse shapes. Fig. 2b then shows the frequency domain plot of the same pulse shapes as 

20 well as the noise power spectral density averaged over 400 simulated noise traces. Both the pulses and the noise 

21 spectral density agree well with that predicted from linear small signal responsivity models (such as that described in 

22 refs. [14,22]) using the same device parameters. 

23 

24 3. Correlated energy-position optimal filter (CEPOF) algorithm 

25 The CEPOF algorithm is an adaptation of the standard digital optimal filter for determining photon energy in single 

26 pixel calorimeters [17]. In a PoST we want to solve for both the energy and the position simultaneously. When a 

27 photon is absorbed at position x y on a one-dimensional absorber, the measured data D can be modelled as a (first 
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order) Taylor expansion around some reference position x such that D = H(S + (. x y - x )S / ), where S is the signal 
responsivity, = dS/dx is derivative of S with x and H is a scalar representing the photon energy. Thus, to construct 
our CEPOF we minimise the difference between the measured data and our modelled pulse. In one-dimension, with 
n TES readout channels, the % 2 condition to minimise is [16]: 


x 1 = X [ D - H i s + &cS')J T w[d~h(s + 8xS‘ )] , 


(i) 


where Z), S are S are n element vectors containing the data, responsivity and derivative of the responsivity 
respectively for each TES. W is the weight matrix and is the inverse of the n x n variance-covariance matrix, which 
describes the detector noise. Eq. 1 is easily adapted for the Hydra case, which consists of a singe TES (n = 1) 
connected to a series of m absorber elements: 


x 2 =E 


(D-H(s,+5iS!jf 

H) 


( 2 ) 


here, TV 2 is the measured power spectral density of the noise. Even though they may not be physically positioned in 
any particular order, for the purpose of this algorithm we can reference the absorbers in order of most similar pulse 
shapes from / = 1 . . . m. Thus we replace the continuous position variable x with an integer position /. Then, 5, is 
calculated for each absorber / and their derivatives are evaluated between sequentially ordered absorbers such thatS^j 
= dSi/di. 8/ is then simply (/ y - i). As described in [16], we can minimise Eq. 2 with respect to 8/ and H , which leads 
to a pair of simultaneous equations in Hhi and //which have the solution: 


H 'l 

[Ha) 




f -X'vs.YXf.o 
l-XYs, 1 f , s , AX Yd 


(3a) 


Collecting normalisation terms into the constants k,\, k l2 and k l3 we can write 
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(3b) 


2 H = k, t J J F l D,+k a Y 4 F;D, 

3 

4 and 

5 

6 H8i = k i2 J j F l D i + k i3 J j F i 'D i . (3c) 

7 

8 For our Hydra we have m template optimal filters of the form F, = S*/N 2 . F/ is then the derivative of the filter with 

9 respect to position z. Thus if we apply a filter corresponding to an absorber z the CEPOF will give a self consistent 

10 estimate of both the energy H and the relative position 8z from the assumed starting position /. If the first order 

1 1 derivatives exactly match the variation in pulse shape across all absorbers then this will give the correct energy and 

12 position. However, in reality the first order terms may only be accurate close to the true absorber position. 

13 Consequently some iteration may be required to converge on the correct pixel. Fig. 3a shows, for the Hydra case 

14 outlined in Section 2, the calculated value of 8/ (using Eq. 3) for each pixel z, for a single 100 eV X-ray with no 

15 noise, absorbed in pixel z = 5. The filters and derivatives are all calculated using the 100 eV template pulse and the 

16 average of 400 noise traces. The presence of noise and error in the derivative function at larger 8 z means that 8/ may 

17 not be an exact integer value. Thus to converge on a solution we must round Sz to the nearest integer value. In this 

1 8 example, the maximum number of iterations after the initial estimate required to converge to 8 z = 0 is 3, though the 

19 typical result is only 1 iteration. Fig. 3b shows the calculated value of the photon energy H for the same test case. 

20 Also shown in this plot is the effect of dropping the correlation term in Eq. 3 F*S (k,2 = 0) then our energy 

21 calculation becomes independent of position and H r E (for z z y ) becomes significantly larger. 

22 

23 4. Energy resolution and position sensitivity 

24 Assuming we have converged to the correct pixel (Sz = 0), the theoretical root-mean-squared (rms) energy resolution 

25 and position resolution are found from the variance on H and 7/8/ [16,23]: 

26 
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(4a) 


and 


A E = 




1/2 

2 [I f aX f /sh 


1/2 


(4b) 


For our segmented absorber example, the units of A / are in pixel number or /. Assuming a Gaussian distribution of 8/ 
around 8/ = 0, we can convert Ai into a probability P, that a photon will be assigned to the correct pixel or its 
neighbouring pixels. The probability of an event with |S/1 > x can be found from the complementary Gaussian error 
function, erfc , 


2 °° 2 

P^xa) = A-^r== J e~ u du = Aerfc(xa ) , 


(5a) 


where, 


16 Cl = 


_1 1 _ 

2V2 A i ' 


(5b) 


For the end pixels (i=0 and / = m) we set/4 = 0.5 and for all other pixels (/ = 1 ...w-1), A = 1 . For x = 1, Eq. 5 gives 
the probability that the photon will be assigned to the incorrect pixel (/ i- i y ), forx = 3 this is the probability that the 
photon will have an assignment error of greater than 1 pixel either side of i r Fig. 4 shows the calculated probability 
of an incorrect pixel assignment (x = 1) using Eq. 5 for all the absorber pixels with a photon energy of £= 50 eV. 
Also shown in this plot is the actual error in a sample of 100 simulated photons in which we calculate the pixel 
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1 assignment using Eq. 3. To demonstrate the intrinsic position error at 6/ = 0 we assume a starting position of / = i y 

2 (circle symbols). Both the theoretically calculated and the actual position error show on average very good 

3 agreement across all pixels. In a real detector where we have no a priori information about the event, we must start 

4 at some fixed pixel, such as the center pixel z = 5 for example and allow the CEPOF algorithm to converge to a 

5 solution. Depending upon the level of noise and the accuracy of the derivative functions this can result in a 

6 convergence error due to a local minimum in 8/. This data is also shown in Fig. 4 (star symbols). The total average 

7 position error is slightly degraded from 33.2% to 34.7%. Fig. 5a shows for pixel 5 only, the theoretically calculated 

8 position error as a function of energy as well as the actual position error for 100 simulated pulses at 20, 50, 100 and 

9 300 eV. In Fig. 5b we show, for the same data at 20, 50 and 1 00 eV, a histogram of the actual pixel assignments. 

10 Thus even if the pixel error is relatively high; the majority of incorrectly assigned events are assigned to the nearest 

1 1 neighbour pixels so the position resolution will not be completely lost. At 50 eV the average number of iterations 

12 (across all pixels) after the initial position estimate is 1.23, this falls to 1.12 at 100 eV and 1.02 at 300 eV. For the 

13 latter there is < 0.1 % position error across all pixels. Depending upon the number of absorber elements the number 

14 of iterations could be reduced by making an initial position estimate using some other measure such as rise-time for 

15 example. 

16 

17 At an energy of 50 eV, it was found that approximately 10 % of the total events, averaged over all absorbers, would 

18 fail to converge on a solution and iterated continuously between 2 pixels where in both cases |Sz] > 0.5. In the 

19 implementation of the CEPOF we set the final pixel assignment to correspond to whichever pixel had the 

20 smallest lS/1 . Of these events approximately 50 % were then correctly assigned. By rejecting these events some 

21 modest improvement in the total position error from 35 % to 33 % can be achieved, at the expense of a 10 % drop in 

22 throughput. 

23 

24 The term F/S is a measure of the intrinsic correlation between energy and position and is required to give a self- 

25 consistent value of H and 5/. In a continuous absorber PoST this energy-position correlation results in a degradation 

26 of the energy resolution and position resolution. In the segmented PoST, however, once we have converged on the 

27 correct pixel, the position-energy correlation disappears and this term can be dropped in calculating the photon 
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1 energy. Eq. 5b then reduces to the single pixel form of the energy resolution given by AE, = (4S/ 7 ,5 , / )"° 5 . Fig. 6 

2 shows the calculated FWHM energy resolution (2.355A£,) using the full CEPOF and that with the correlation term 

3 dropped. Thus dropping this term gives a more accurate estimate of the energy. At lower photon energies where we 

4 may have error in the position, then the resolution calculated using the full algorithm gives a better approximation to 

5 the actual energy resolution. 

6 

7 So far the analysis has assumed that the template responsivities and derivatives are calibrated in the small signal 

8 limit where the pulses are linear with energy. Assuming a higher calibration energy such as that from an 55 Fe source 

9 of 5.9 keV has the potential to introduce some level of non-linearity in the detector response. This has the potential 

10 to degrade the position sensitivity since the high energy template may not be well matched to the low energy linear 

1 1 pulse. The vast majority of non-linearity comes about from the dis-proportionality between the measured change in 

12 current A I and the TES resistance R (AI 1 /R 2 ); this effect is intrinsically included in our model. The significance of 

13 this effect will depend upon the specific design parameters of the device in question. For our current Hydra example 

14 we find that by using the non-linear 5.9 keV template, we introduce no additional error in the position sensitivity 

15 though a small correction is required to re-calibrate the gain scale. At higher energies further non-linearity can occur 

16 due to the dependence of a and [5 on the resistive transition and saturation can occur when the TES is driven fully 

17 into the normal state. This effect does represent an extreme case and is neglected in our model since we keep a and 

1 8 fi constant. If non-linearity effects are found to be a problem for a specific detector then using multiple calibration 

19 lines, or implementing non-linear optimal filtering techniques as discussed in ref [24] can be used to improve 

20 performance. 

21 

22 PoST’s have an intrinsic trade-off between energy resolution and position sensitivity, which depends upon the ratio 

23 of the internal to external thermal conductances [16]. If the absorbers are all strongly thermally coupled to the TES, 

24 then the energy resolution will be optimal but there will be no signal variation with position. Reducing the internal 

25 conductances relative to the bath conductance introduces position dependence but attenuates the decoupled absorber 

26 signal relative to the noise (see Fig. 2), as well as adding more internal thermal noise. This is illustrated in Figs. 7 

27 and 8 were we show the calculated FWHM energy resolution and position error respectively at 50 eV as a function 
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1 of the internal link conductances, which are normalised to the configuration described in Table 1. Only a subset of 

2 pixel data are shown. As the link conductances are reduced the position error passes through a minimum and then 

3 starts to degrade. This is because the degradation in signal-to-noise out-weighs any improvements due to increasing 

4 the separation of the pulse shapes. The minimum occurs at different conducatances for each pixel but are relatively 

5 broad, meaning it is possible to achieve both good position sensitivity and energy resolution across all absorber 

6 pixels. As the G//s are reduced, the TES and most strongly coupled absorber become more and more isolated from 

7 the other elements and the energy resolution (without position correlation) improves. In the extreme limit of Gi -» 0 

8 we find A£ = 1.7 eV; since the energy resolution scales as the square root of the total heat capacity this is V 9 times 

9 lower than the G L — > °° case where A E =5.1 eV. 

10 

1 1 5. Effect of random arrival phase 

12 A potential difficulty with the implementation of the digital optimal filter on real data is that the measured signal is 

13 random in time and its correct phase relative to the filter is unknown. The variation in the arrival phase will depend 

14 upon the noise level relative to the signal and the shape of the initial pulse rise, as well as the method of triggering 

15 and the sample rate. Although calculated in the frequency domain, the optimal filter is usually applied as a cross- 

16 correlation in the time domain. If the phase of the filter and the pulse are the same then this will correctly give the 

17 photon energy (as well as position for the CEPOF). If the arrival phase is not known then the complete convolution 

18 over all phases is required. In the simple OF the peak in the convolution gives the correct arrival phase and photon 

19 energy with the lowest noise. Close to the maximum, phase errors give small errors in energy but if the phase error 

20 is large then the error in energy can become very much greater than the actual noise. Thus for discreetly sampled 

21 data it is typical to apply some phase correction algorithm, such as interpolating the energy between the quantised 

22 phase lags around the peak to more accurately measure the energy. Fig. 9a shows a plot of the convolution response 

23 as a function of phase lags for a 300 eV (noise free) pixel 1 pulse, calculated using both the CEPOF and the standard 

24 OF. Unlike the standard OF, the CEPOF convolution is asymmetric about zero lags and the peak cross-correlation 

25 does not give the optimum noise or the correct energy. This is the same as noted by Whitford (2005) [24] for similar 

26 non-linear (in energy) optimal filtering. The standard OF is sharply peaked around 0-lags because of the high 

27 frequency content in the filter and the pulse, and relatively small errors in lag could create a large error in energy. 

28 Fig. 9b shows the case for a pixel 8 pulse. For this slower pulse the dependence on lags is much less, so the potential 
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1 for error is smaller. Fig. 9c then shows 8 i for pixel 1 and 8 using the CEPOF. Since we only require | 8*1 < 0.5 for 

2 the correct pixel assignment, small phase shifts around the true phase do not result in any position error. If the jitter 

3 is greater than a few lags then some correction will be required. The fact that the peak in the convolution for the 

4 CEPOF does not correspond to the correct energy is an obvious problem for the implementation of the algorithm on 

5 random discretely sampled data. However, it is possible to correct for phase (at the level of integer lags) using the 

6 standard OF first and then, by applying the CEPOF we can converge on the final position. A final sub-lag phase 

7 correction (by interpolating between lags) using the standard OF can be used for the final estimation of the energy. 

8 

9 6. Summary and conclusions 

10 We have demonstrated an implementation scheme for a correlated-energy-position optimal filter algorithm, which 

1 1 we have adapted for single-TES, segmented-absorber Hydras. Although in this paper we demonstrate the algorithm 

12 on Hydras devices, we have also implemented the algorithm on simulated 2-TES segmented PoST data with similar 

13 results. The CEPOF can be extended to include extra dimensions and any number of channels for segmented or 

14 continuous absorber devices. The CEPOF can also be implemented on other calorimeter technologies such as 

15 position-sensitive MMC’s [18]. For a continuous device the CEPOF implementation will in principle be the same, 

16 however we have a continuous distribution of pulse shapes rather than a quantised set. The requirements on 

17 calibration are potentially higher for continuous devices, though could be reduced by interpolating the filter 

1 8 functions between a set of pre-computed calibration points at known xlL. The computational requirements of the 

19 CEPOF are only moderately more than the standard OF. Every absorber element whether in a single pixel array or a 

20 PoST array requires its own template filter. For the CEPOF we require the calculation of the additional 

21 normalisation constants (which only need to be calculated once) and the extra derivative filter convolution term. On 

22 the simulations reported here we on average required only one additional iteration of the CEPOF to converge on the 

23 correct position, hence the CEPOF is not likely to be significantly more computationally intense than the standard 

24 OF. 

25 
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Table 1. 9-Pixel Hydra model parameters. Note that absorber one is assumed infinitely well coupled to the TES. 


Parameter 

Symbol and units 

Value 

Heat capacity per absorber 

C (pJ/K) 

0.35 

Thermal conductance to heat sink 

G h (pW/K) 

650 

TES transition temperature 

T e (mK) 

90 

Bath temperature 

T b (mK) 

50 

TES normal state resistance 

R-n (mQ) 

7 

Fractional bias point 

R 0 /Rm (%) 

15 

Bias circuit load resistor 

Ri (m£2) 

0.25 

Internal thermal conductances for absorbers 2-9 

G l (nW/K) 

34,21, 14, 10.3, 
7.3, 5, 3.7, 3 

TES resistance sensitivity 

a 

60 

TES current sensitivity 

p 

2.0 

SQUID input coil inductance 

U n (nH) 

90 

Integration step size 

At (ps) 

1.0 

Number of samples per record 

N, 

T* 


15 



1 Figure captions 

2 Fig. 1. (color online) Schematic diagram of the 9-pixel Hydra. 

3 Fig. 2. a) (color online) Simulated 9-pixel Hydra noise-less pulse shapes for a photon energy of 100 eV. b) Power 

4 spectral density averaged from 400 simulated noise records and spectral density of the noise free pulse shapes. 

5 Fig. 3. a) (color online) Calculated value for 8/ at each of the 9 Hydra absorber pixels for a 100 eV photon absorber 

6 at pixel i = 5. b) For the same example we show the calculated photon energy H using the complete algorithm and 

7 calculated assuming position and energy are independent. 

8 Fig. 4. (color online) Theoretically calculated position error for each pixel at a photon energy of 50 eV. Also shown 

9 is the actual position error for a population of 100 simulated pulses assuming a starting position of / = / x and then for 

1 0 the case / = 5 . 

1 1 Fig. 5. a) (color online) Theoretically calculated position error as a function of energy for pixel 5. Also shown is the 

12 actual position error for a population of 100 simulated pulses at energies of 20, 50, 100, 300 eV. b) Histogram of 

13 pixel assignments for simulated data at 20, 50 and 100 eV. 

14 Fig. 6. (color online) FWHM energy resolution calculated using both the full CEPOF algorithm and dropping the 

15 energy-position correlation. 

16 Fig. 7. (color online) FWHM energy resolution calculated using both the full CEPOF algorithm (solid lines) and 

17 dropping the energy-position correlation (dashed lines) as a function of the internal link thermal conductances G L 

18 ( G l is normalized to configuration described in Table 1. Data is shown for pixels 1, 5 and 8. 

19 Fig. 8. (color online) Position error at 50 eV as a function of normalized link conductance for pixels 1, 2, 5 and 8. 

20 Fig. 9. a) (color online) Convolution of the CEPOF and standard OF with a 300eV pulse from pixel 1 as a function 

21 of phase lags. Zero lags corresponds to the correct arrival phase of the pulse relative to the filter, b) Same as for Fig. 

22 9a using a pixel 8 pulse, c) 8/ for both pixel 1 and pixel 8. 
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